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ABSTRACT 

We present the results of our study of spectral energy distributions (SEDs) of a sam- 
ple of ten low- to intermediate-synchrotron-peaked blazars. We investigate some of 
the physical parameters most likely responsible for the observed short-term variations 
in blazars. To do so, we focus on the study of changes in the SEDs of blazars corre- 
sponding to changes in their respective optical fluxes. We model the observed spectra 
of blazars from radio to optical frequencies using a synchrotron model that entails a 
log-parabolic distribution of electron energies. A significant correlation among the two 
fitted spectral parameters (a, b) of log-parabolic curves and a negative trend among 
the peak frequency and spectral curvature parameter, b, emphasize that the SEDs of 
blazars are fitted well by log-parabolic curves. On considering each model parameter 
that could be responsible for changes in the observed SEDs of these blazars, we find 
that changes in the jet Doppler factors are most important. 
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1 INTRODUCTION 

Studies of the two types of blazars, BL Lac objects and 
flat spectrum radio quasars (FSRQs), have shown that their 
Spectral Energy Distributions (SEDs) are characterized by 
a dou ble peaked luminosity structure (e.g. iGhisellini et al.l 
1 19971 ). The BL Lacs and FRSQs classes are deflned ac- 
cording to the absence or presence of strong broad emis- 
sion lines in their optical/UV spectrum, respectively. The 
shapes of these bumps are characterized in the log{i/F^) 
vs. log v plot by a smooth spectrum extending through 
several frequency decades. The low-energy peak is well ex- 
plained by synchrotron emission from relativistic electrons 
in a jet closely aligned to the line of si ght, with bulk 



Lorentz factor, F of order of 10-100 



e^. 



Maraschi et al.l 



Il992l : IGhisellini et al.lll993l : iHovatta et aDl2009l . and refer 
ences therein). Moreover, the position of the low-energy peak 
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leads to a further classification of the blazars into three cate- 
gories, depending on the peak frequency of their synchrotron 
bump, Vp^ak'- low-synchrotron-peaked (LSP) sources with 
^peafe < 10^^ Hz; intermediate-synchrotron-peaked (ISP) the 
sources, those with 10^* Hz < v^^ak < 10^^ H^i ^'^id high- 
synchrotron-peaked (HSP) sources, with Vp^a,k > 10^^ Hz. 
T his scheme is an extension o f the classificatio n introduced 
by lPadovani &: Giommil (| 19951 ) for BL Lacs (see lAbdo et al.l 
i2010bl . for details'). 

The high energy (hard X-ray and 7-ray component of 
the SED is usually explained as arising from inverse Comp- 
ton scattering of the same electrons producing the syn- 
chrotron emission. These electrons interact either with the 
synchrotron photons th emselves (synchrotron self-Compton, 
SSC; e.g., Marsch er fc G ear 1985) or with external pho- 
tons originating in the local environment (external Comp- 
ton, EC). In the latter case soft ph otons can be produced 
direc tly by the accretion disk (e.g. iDermer fc Schlickeisen 
[1993), or indirectly, for instance those reprocessed by the 
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broad line region (BLR) (e.g. ISikoralll994l '). or by the dust 
torus (|Blazeiowski et al.ll2000l ). Alternative hadronic mod- 
els, where the 7-rays are produced by high-energy protons, 
either via proton synchrotron radiation or via secondary 
emission from photo-pion and ph oto-pair-production reac- 
tions, have also been proposed (see lBottcher et aUbOOTl . and 
the references therein for a review). 

The shape of these bumps is characterised in the 
hogiyFv) vs. Log v plot by a smooth spectrum extend- 
ing through several frequency decades. Below and above 
the peaks the spectrum can be approximated with simple 
power l aw profiles with spect ral indices above and below the 
peaks (|Ghisellini et al.lll99a ). and such power law spectra 
are naturally produced if the the emitting electrons follow a 
power law distribution of energy. So, the broad band SED 
of blazars can be well approximated by a simple parabolic 
funct i on wi t h logarithms of its vari ables (e.g.,|Massaro ct al. 
|2004 I2OO6I : iTramacere et al.ll2007l . and references therein). 
The log-parabolic function is one of the simplest ways to 
represent curved spectra and under simple approximations, 
can be obtained via a statistical electron acceleration mech- 
anism where th e acceleration probab ility decreases with the 
particle energy (iMassaro et al.ll2006l ). 

The emission from blazars is known to be variable at 
all wavelengths. The flux variability is often accompanied 
by spectral changes. These SED changes are very likely as- 
sociated with changes in the spectra of emitting electrons 
that arise from differences in the physical parameters of the 
jet. Hence modelling of blazar broadband spectra is required 
to understand the extreme conditions within the emission 
region. Not only is the broadband SED crucial to this un- 
derstanding, but variability information is also needed to 
allow us to describe how high emission states arise and how 
they difi^er from the low states. This type of study is most 
important in discriminating between models. Since it is rea- 
sonably to assume that only a few parameters of the model 
will change significantly between different states of the same 
source, such comparative modelling of broadband spectra 
allows us to put rather tight constraints on those model 
parameters that are likely to change, at least under the as- 
sumption that all other parameters are held flxed for the 
different model fits (e.g. Mukheriee et al.lll999l : IPetrv et al.l 



I2OOOI : iHartman et al.|[200ll i 



In this paper, we focus on the study of changes in the 
SEDs of BL Lacs and FSRQs corresponding to observed 
changes in their respective optical fluxes. The SED varia- 
tions are expected to be produced by changes in the spectra 
of the emitting electrons which in turn arise from variations 
in the physical parameters of the emission region. There- 
fore we needed to first construct models and then inves- 
tigate how the SEDs observed in different states might be 
explained through changes in these parameters. More explic- 
itly, we modelled the observed spectra of blazars from radio 
to optical frequencies using a synchrotron model with the 
emitting electrons following a log-parabolic distribution of 
energy. This allows us to try to estimate the factors respon- 
sible for short term optical flares or short-term variability 
(STV) in blazars. In this paper, we will strictly concentrate 
on the low energy hump (synchrotron emission) portion of 
blazar SEDs. A complete investigation of the entire broad 
band SEDs of these objects, including X-ray and 7-ray data, 
will be performed in the future. 



This paper is structured as follows. Section 2 provides 
a brief description of the multi-frequency data we employed. 
In Section 3, we discuss the SED modelling and Section 4 
provides our results. Our discussion is given in Section 5 and 
we present our conclusions in Section 6. 



2 MULTI-FREQUENCY DATA 

The multi-frequency SED data of our 10 blazars span a fre- 
quency range between radio and optical, including mm, sub- 
mm and infra-red (IR) data. These measurements were all 
taken between 2008 September and 2009 June. Our sample 
consists of five BL Lacs and five FSRQs, eight of which are 
LSPs and two of which are ISPs. We selected these sources 
as they were bright enough for us to obtain day to day flux 
coverage in B, V, R, and I passbands. Certainly this num- 
ber of objects is too small to provide clearly representative 
samples of the BL Lac and FSRQ classes and there may be 
unknown biases in our two groups. Despite the small sam- 
ple sizes, some general behaviours might be noted in these 
data and then investigated in more detail in the future with 
larger samples. The short term variability (STV) in flux as 
well as in colour o f these blazars was recently reported by 
I Rani et al.l l|2010bl ). We used this optical data in our SED 
study and added to them data collected at other frequencies 
over the same time period. Radio fiux densities at 4.5 GHz, 8 
GHz and 14.5 GHz frequencies are obtained from University 
of Michigan Radio Astronomy Observatory (UMRAQj) data 
base. The mm and sub-mm data are provided by the Sub- 
millimeter Array (SMA) Observer CenteiQ data base. The 
near-IR data are collected from the monitoring provided by 
the Small and Moderate Aperture Research Telescope Sys- 
tem (SMART© . 

For investigating the impact of the variations in flux 
on the source parameters used for modelling the SEDs we 
always ascertained two SEDs of the same source that are 
characterized by a change of at least 0.3 magnitudes in the 
optical R-band. The brighter one is denoted as the high-state 
and the fainter as the low-state. The time periods during 
which the different SEDs of all the sources were obtained 
are listed in the first column of Table 1 immediately be- 
low the name of each of t he blazars. Details of their optical 
properties can be found in I Rani et al.l (|2010a ). We note that 
while we tried to obtain and utilize data taken simultane- 
ously, this was rarely possible, and the temporal intervals 
over which the data needed to produce an SED in a low or 
high state could be found range from three days through 
three months. Clearly, this lack of simultaneity can lead to 
substantial uncertainties in the analysis. 



3 SED MODELLING 

We performed the spectral analysis using a homogeneous 
synchrotron emission model with log parabola (LP) energy 
distribution of emitting electrons to fit the lower energy part 
of observed spectra, i.e., synchrotron spectra. The best fit 



^ http://www.astro.lsa.umich.edu/obs/radiotel/umrao.php 

^ http://smal.sma.hawaii.edu/callist/callist.html 

^ http://www.astro.yale.edu/smarts/glast/targets.html 
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model to each blazar in each of two states was obt ained by a 
numerical SSC code (JTramacere et al.ll2007l . [20091 ). We used 
an SED codqj available on-line. 

The simplified model assumes that radiation is pro- 
duced within a single blob in the jet, which is taken to be 
moving relativistically at small angle along the line of sight 
of the observer. Thus the observed radiation has a Doppler 
boosting factor 5 = [r(l — j3 cos8)]~^ , where /3 is the ve- 
locity of the source divided by the velocity of light, F is the 
Lorentz factor, and 6 is the angle between the line of sight 
of observer and direction of motion of the source. The ob- 
served SED of blazars has a double peaked structure and a 
simple analytical function that can model the shape of these 
broad peaks is a parabola in the logarithms of the variables, 
i.e., a log parabola (LP). This fu nction has three spe ctral 
parameters and can be defined as (jMassaro et al.ll2006l ') 



F{E) = KE 



-{a+b log{E)) 



keV 



(1) 



where a is the photon index at an arbitrary energy (usu- 
ally taken to be 1 KeV) and b is the measure of spectral 
curvature of the observed radiation. One advantage of the 
LP functional form compared to other option is that in LP 
formulation the curvature around the peak is characterized 
by a single parameter b, while in other models it is char- 
acterized by more co mplex functions (e.g. Sohn et al. 2003; 
iMassaro et al.l 120061 ): however, the limitation is that it is 
more suited for distributions symmetrically decreasing with 
respect to peak frequency, than is the norm for synchrotron 
spectra. Although more complex models, such as a power- 
law plus log-parabola model, should fit the mm and longer 
wavelengths data better, we do not have good coverage in 
those bands; further, the simpler log-parabola model has a 
clea r physical moti v ation, which we now discuss. 

iMassaro et al.l (|200J) showed that the log-parabolic 
spectrum is analytically related to the statistical acceler- 
ation mechanism. It is obtained when the acceleration prob- 
ability, p, depends upon energy itself, i.e.. 



p^^gh! (i = o,i,2,...). 



(2) 



where g, and q are positive constants. Such a situation can 
naturally occur when particles are confined by a magnetic 
field that has a confinement efficiency that decreases with 
an increase in t he gyration radii of the accelerating particles 
iJMassaro et al.ll2006l ). The integrated energy distribution of 
such accelerated particles is given by 

r log(-t/'iQ)] 



N{> 7) = iV(7/7o)"'' 



(3) 



where 70 is the minimum Lorentz factor in Eq. (2), and r 
and s are the spectral parameters of the electron popula- 
tion and are related to those of the emitted radiation as 
IMassaro et al.ll2006h 



l)/2 



and 



r/4. 



(4) 



With these approximations, we can completely specify 
the SSC model with the following parameters: magnetic field 
intensity (B), size of emission region (R), Doppler boosting 
factor (5), the log-parabola spectral indices (s and r), the 
number density of emitting electrons (N), the redshift of the 
source, z, and particle Lorentz factors, "/min, 'y-max and 70. 

* http://tools.asdc.asi.it/ 



While fitting the SED of blazars using these parameters, 
we know the value of z and we hold R, "fmin, and jmax 
fixed. It is clear that there is a cut-off in the low energy part 
of synchrotron spectra at a frequenc y around 10^^ Hz due 
to synchrotron self-a bsorption (e.g. iGhisellini et al.l 1 1993 : 
iTavecchio et al.l |2002| . and references therein), so we con- 
sider this self-absorption in modelling the SED of blazars. 
Because we are looking at short term variability over the 
course of a just a few months at most, we feel it justified 
to assume that R is essentially constant for each source and 
we find that the values of ymin, and ^max we used for every 
case give good ranges for the electron energies needed to fit 
any of the SEDs and so they are fixed as well. The values of 
r, s and 70 are rather tightly constrained to fit the observed 
slopes and we then made over than 30 different models for 
each source to produce the values of all the parameters that 
yield a "best fit". 

A direct observation of the synchrotron peak will give a 
significant help to constrain the model parameters, but un- 
fortunately, observations of blazars in infrared bands, above 
all, in the medium and far infrared, are not so common as 
are the radio, optical and near IR ones. Occasionally this gap 
can be covered by s pace-based telescopes such as Herschel 
(|Pilbratt et al.ll2010l i. In addition, for FSRQs the SEDs are 
sometimes overwhelmed by thermal emis sion from an accre- 
tion disk, particularly during low-states (jMalkan fc Sargenj 
11982 !'), but our code is not able to consider this contribution 
while modelling the SEDs of blazars. This additional contri- 
bution may explain the excess of optical/UV emission with 
respect to the best-fit model for some FSRQs. The values 
of all these parameters that provide good fits for the two 
different SEDs of all of the eight LSPs and two ISPs in our 
sample are listed in Table 1. The fitted SEDs are shown in 
Figs. 1 - 3. 



4 RESULTS 

We fit the radio to optical SEDs of 10 blazars with a syn- 
chrotron model. The SED curves of all the sources are shown 
in Fig. 1 — 3 and the fitted parameters are listed in Table 
1. Detailed multiband optical STV studies of the fluxes and 
colours of a ll of these blazars o ver the same time period is 
reported in lRani et al.l (|2010bl ). There we showed that the 
colour versus brightness correlations support the hypothe- 
sis that these BL Lacs tend to be bluer with an increase in 
brightness, while these FSRQs show the opposite trend. Our 
intraday variability (IDV) study of these sources in the opt i- 
cal R-band have been recently reported in lRani et al.l l|2011f ). 
All of these blaz ars belong to the F irst Fermi LAT catalogue 
(IFGL Catalog. I Abdo et al.ll2010bl ). and their spectral prop- 
erties (hardness, curvature and variabil i ty) at GeV energies, 
have been established bv lAbdo et al.l (|2010d V The distri- 
bution of photon indices (F) above 100 MeV is found to 
correlate strongly with blazar subclass. Also, the spectral 
indices tend to be harder when the flux is brighter for both 
FSRQs and BL Lacs. 

We now summarize some previous observations of 
each of these sources. Since the optical flux, colour and 
spectral variabilit y of these LSPs a nd ISPs have already 
Rani et al.] (|2010bf ). and IDV studies are 



be en discussed in 

in I Rani et al.l (|201ll ). here we highlight some results from 
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Table 1. Parameters of SED models using log-parabola fits 



Source 


< Rmag > 


Log R 


B 


S 


N 


z 


Log 


Log 


Log 


r 


s 


Log vF^pcak 


Log vpeak 






(cm) 


(Gauss) 




(cm-3) 




Imin 


Imai 


70 






(org cm^^ s^^) 


(Hz) 


3C 66A 




























20-30 Oct. 2008 


14.60 


17.0 


0.10 


15.2 


30 


0.444 


0.2 


6 


4.00 


0.5 


3.00 


-10.71 


14.66 


20-22 Jan. 2009 


14.05 


16.8 


0.11 


12.2 


30 


0.444 


0.2 


6 


4.10 


0.6 


3.00 


-10.40 


14.80 


AO 0235-1-164 




























Sept. 2008 


16.00 


17.2 


0.09 


20.0 


40 


0.94 


0.2 


6 


3.38 


1.9 


3.65 


-10.21 


13.07 


20-30 Oct. 2008 


15.00 


17.2 


0.09 


21.0 


40 


0.94 


0.2 


6 


3.42 


1.9 


3.65 


-10.07 


12.95 


20-22 Jan. 2009 


17.00 


17.2 


0.08 


18.0 


35 


0.94 


0.2 


6 


3.34 


1.9 


3.40 


-10.58 


13.03 


PKS 0420-014 




























20-30 Oct. 2008 


16.85 


17.0 


0.60 


18.0 


35 


0.915 


0.2 


6 


2.80 


1.1 


3.10 


-10.81 


12.93 


Jan. 2009 


16.45 


17.0 


0.60 


20.0 


37 


0.915 


0.2 


6 


2.80 


1.1 


3.20 


-10.69 


12.95 


S5 0716-1-714 




























20-30 Oct. 2008 


15.30 


17.0 


0.09 


14.0 


20 


0.31 


0.2 


6 


3.60 


0.9 


3.20 


-10.78 


13.60 


Mar.-May 2009 


13.60 


17.0 


0.085 


15.0 


15 


0.31 


0.2 


6 


4.12 


0.75 


3.34 


-10.29 


14.55 


PKS 0735+178 




























May 2009 


15.90 


17.13 


0.07 


12.0 


20 


0.424 


0.2 


6 


4.00 


0.5 


3.00 


-11.20 


14.41 


June 2009 


15.60 


17.13 


0.07 


12.7 


20 


0.424 


0.2 


6 


4.00 


0.5 


3.02 


-11.14 


14.52 


OJ 287 




























20-30 Oct. 2008 


14.20 


17.0 


0.10 


18.0 


20 


0.306 


0.2 


6 


3.40 


1.5 


3.25 


-10.07 


13.51 


Mar.-Apr. 2009 


14.70 


17.0 


0.12 


19.0 


25 


0.306 


0.2 


6 


3.30 


1.5 


3.45 


-10.06 


13.21 


4C 29.45 




























20-29 Mar. 2009 


17.00 


17.0 


0.50 


15.0 


30 


0.729 


0.2 


6 


2.80 


1.3 


2.95 


-10.77 


12.96 


May 2009 


18.20 


17.0 


0.45 


15.0 


25 


0.729 


0.2 


6 


2.80 


1.1 


2.95 


-11.16 


12.91 


3C 279 




























Apr. 2009 


15.40 


17.0 


0.60 


17.0 


30 


0.5362 


0.2 


6 


2.70 


1.7 


3.20 


-10.20 


12.82 


June 2009 


16.00 


17.0 


0.65 


17.0 


30 


0.5362 


0.2 


6 


2.70 


1.4 


3.30 


-10.43 


12.53 


PKS 1510-089 




























Apr. 2009 


15.40 


16.8 


0.50 


14.0 


30 


0.36 


0.2 


6 


3.10 


1.0 


3.155 


-10.69 


13.41 


June 2009 


16.00 


16.8 


0.50 


14.0 


30 


0.36 


0.2 


6 


3.10 


0.8 


3.155 


-11.02 


13.45 


3C 464.3 




























4-10 Sept. 2008 


14.70 


17.2 


0.50 


23.0 


30 


0.859 


0.2 


6 


2.80 


2.3 


3.35 


-9.46 


13.22 


20-30 Oct. 2008 


15.40 


17.2 


0.48 


21.0 


30 


0.859 


0.2 


6 


2.70 


1.9 


3.00 


-9.81 


12.75 



Log R(cm) : Size of emitting region 

B(Gauss) : Magnetic field 

3 : Doppler boosting factor 

N(cm~ ) : Number density of emitting electrons 

z : redsliift 

"fmim "fmax '■ Minimum and maximum values of Lorentz factor 

7Q : tbe Lorentz factor corresponding to the energy where s is evaluated 

r and s : spectral parameters of the electron population (see text) 

^Peak '■ Synchrotron peak frequency in the rest frame of source 



earlier studies of these sources at other wavebands. 

3C 66A 

This ISP blazar is classified as a BL Lac object. The redshif t 

of 3C 66A was reported as z=0.444 

but this value remains uncertain (see iBramel et al. 



197g) 



l2005l : 



lYan et al.ll2010l ). 3C 66A is observed in radio, IR, optical 
X-rays, an d Kamma rays and sh o ws strong luminosit; 



variations I Terasranta et al.1 |2004| : iBottcher et al. 200' 



I Rani et alj l2010bl . and references therein) . iBottcher et al.l 



I 



1 2005^ organized an intensive multiwavelength campaign of 
the source from 2003 July through 2004 April to understand 
its broadband spectral behaviour. Ijoshi fc Bottcheij (120071 ) 
used a leptonic SSC jet model to reproduce the broadband 
SED and the observed optical spectral variability patterns 
of 3C 66A and made predictions regarding observable X-ray 
spectral variability patterns and 7-ray emission. 

AO 0235+164 

The blazar AO 02354-164 at z = 0.94 iJNilsson et all 119961 ) 
was classified as a BL Lac object by Spinrad & Smith 
(1975). An apparent quasi-periodic oscillation of '-^17 days 
has been repor te d in the X-ray Ught curves of the source 
by iRani et al.1 (|2009l ). The XMM- Newton observations 
suggested the presence of an extra emission component 
in the source SED, in addition to the synchrotron and 



inverse-Compton ones. The origin of this component, 
peaking in the UV/soft X-ray frequency range, is probably 
thermal emission from an accret ion disc, or synch rotron 
emission from an inner jet region (IRaiteri et al. 2006, |2008| . 
and references therein) . iHagen-Thorn et al.l (|2008i ) reported 
a significant correlation between the flux density and degree 
of polarization at optical frequencies and found that at the 
maximum degree of polarization the electric vector tends 
to align with the parsec-scale radio jet's direction. 

PKS 0420-014 

PKS 0420-0 1 4 is a FSRQ at a redshift z = 0.915 
IjKuehr et al.l Il98ll ) that has been observed in optical 
bands since 1969. There is a series of papers reporting the 
optically active and bright phases of source and semi-regular 
major flaring cycles (e .g. IWebb et al.l Il988l : IVillata et al 



1997 1 Raiteri et al.lll998l , and references therein) . iRani et al.. 
( 2010br ) reported that the source tend to be redder with 



increase in its brightness. Recently rapid optica l and radio 
bright ening of the blazar has been reported by iBach et al.l 
(|2010h . 



S5 0716+714 

S5 0716+714 was id entified as an ISP JGiommi et al.]ll999l ). 
at z = 0.31 + 0.08 iJNilsson et aIll20od ). However. Nieppola 
et al. (2006) studied the SED distribution of a large sample 
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Figure 1. The modelled SED curves of the blazars 3C 66A, AO 0235+164, PKS 0420-014 and S5 0716+716. The points represent the 
observed data while the best fitted model is shown by line. Since the low frequency part (:?; lO'^^ Hz) of blazars SEDs is governed by 
synchrotron-self absorption mechanism, so the modelled SED is steeper below a frequency of 10^^ Hz (see text for details). 



of BL Lac objects and categorized 0716+714 as a LSP. Two 
strong 7-ray flares on September and Octo ber 2007 from 
the source have been detected by AGILE (|Giommi et al.l 
l2008l : IChen et al.ll2008h . They have also carried out the SED 
modelling of source with two synchrotron self-Compton 
(SSC) emission models, representative of a slowly and 
a rap idly variable component, respectively. IVillata et al.l 
l|2008l ) reported optical-IR SED study of the s ource during 
the G ASP-WEBT-AGILE campaign in 2007. iGupta et all 
l|2009l ) found high probabilities of quasi-periodic com- 
ponents of ~ 25— 73 minutes in the optical light curves 
observed by iMontagni et al.l (|200y). Recently, nearly 
periodic oscillations of ^15 minutes have been r eporte d 
in the optical light curve of the source bv lRani et alj (|2010al ) . 



PKS 0735+178 

PKS 0735+718 is a highly variable q uasar and belong s 
to the category of BL Lac objects (Carswell et alJ 11974 '). 
There have been s everal papers conce r ning i ts rcdshift 
deter mination (e.g. ICarswell et al.l Il974l : iFalomo fc Ulrichl 
|2000| . and the references therein) with the most recent 
resuh of z = 0.424 for PKS 0735+718 found using a Hubble 
Space Telescope (HST) snapshot image (Sbarufatt i et al.l 
[2OO5I); using that redshift, we classify this source as an ISP. 
The synchrotron part of SED of the source peaks at IR-UV, 
while the low X-ray variability with respect to the high 
optical-IR variations, supports the idea that X-rays are pro- 
duced by an inverse Compt on mechanism (Bre gman et al.l 
1 1984 iMadeiski fc Schwartd Il988l ). Using ten yea r loni 
multiband optical data of this source ICiprini et al.l (|200 
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Figure 2. The modelled SED curves of the blazars PKS 0735+178, OJ 287, 4C 29.45 and 3C 279. 



studied the multicolour behaviour, spectral index variations 
and their correlations with brightness. 

OJ 287 

OJ 287 (0851-+202) (z = 0.306) is among the most ex- 
tensively observed and best studied BL Lac objects with 
respect to variabihty. As it is very bright, it is among the 
very few AGN's whose optical observations are available 
for m ore than a century (|Sillanpaa et al .11 199 a : I Gupta et al.l 
I2OO3). A binary black hole model has be en proposed to 
expla i n the ~12 years opt ical periodicity (jSillanpaa et al.l 
119961 : IValtonen et al.l 120081 ). An intense optical, infrared, 
and radio variability study of the source for a time period 
betwe en 1993 to 1998 was been reported in iPursimo et al.l 
IJ2OOC1 I. They detected two major optical outbursts in 
November 1994 and December 1995 even though the radio 
flux was very low during the period of outbursts. Using 



the multi-waveband flux and linear polarization observa- 
tions and combini ng with sub-mm polarimetric images, 
lAgudo et al.l (l201l|) argued that the location of the 7-ray 
emission in prominent flares is > 14 pc from the central 
black hole. 

4C 29.45 

4C 29.45 (z = 0.729; IWills et al.l Il983h . belongs to the 
category of FSRQs. Short and long term variability at 
both IR and optical b a nds has been obser ved in this source 
l|Noble fc Milled Il996l : I Ghosh et al.ll2000l ). Its o ptical flux 
and multi-colour variations were studied by iFan et al.l 
(2006). They reported amplitude variations of ~4.5— 6 
mag in all passbands (U, B, V, R, I) and also suggested 
that there were possible periods of 3.55 and 1.58 yr in 
the long-term optical LC of the source. The short term 
multi-band optical colour-brightness studies showed that 



Spectral Energy Distribution variation in BL Lacs and FSRQs 7 



the source tend to be redder with increase in its brightness 
(jRani et al.1 l2010a ). V ery recently, a GeV flare has been 
detected in this blazar (|Ciprinill20ld ). 



3C 279 

3C 279 is both one of the most extensively observed and 
most violent variable FSRQs across the whole electromag- 
netic spectrum; it is at a redshift oi z = 0.536. A comparison 
of radio to Gamma-ray SEDs of the source during low 
and high states of activity in the WEB T campaign of 
2006 is presented in iBottcher et al.1 l|2007l 'l. ICollmar et al.1 
l|2010l) carried out multifrequency variability study of 
source during January 2006. A complete compilation of 
all simultaneous SEDs of 3C 279 collected during the 
lifetime of CGRO, including modell ing of those SEDs wit h 
a leptonic jet model, i s presen ted in lHartman et al.l (l2001f ). 
Recently, lAbdo et al.l (|2010al 'l report the coincidence of a 
7-ray flare with a dramatic change of optical polarization 
angle. This provides evidence for co-spatial of optical and 
7-ray emission regions and indicates a highly ordered jet 
magnetic field. 



PKS 1510-089 

PKS 1510 -089 is classified as a FSRQ at a redshift of 
z = 0.361 (|Thompson et al.lll99d ). It also belongs to the 
category of highly polarized quasars and already detected 
in the MeV-GeV energy band by the EGRET instru- 
ment on board the Compton Gamma-Ray Observatory 
l|Hartman et al.l Il992l '). The inverse Compton component 
dominated by the 7-ray emission, and the synchrotron 
emission peaked around IR frequencies, even if there is 
clearly visible in this source a pronounced UV bump, 
possib ly caused by the therma l emission from t he ac cretion 
disk (JMalkan fc Moorel Il986l : iPian fc Tre^ Il993l ). The 
multi-wavelength flux and SED variability study of the 
source during its high 7-ray activity period b etween 2008 
Septem ber and 200 9 June has been reported by I Abdo et al.l 
l|2010al 'l. Recently iD'Ammando et al.1 HoTJ) reported the 
detailed analysis of multi-frequency coverage of extreme 
gamma-ray activity from the source observed by AGILE 
in March 2009. They have also carried out broad band 
SED study of the source and found thermal features in the 
optical/UV spectrum of the source during high 7-ray state. 



3C 454.3 

30 454.3 is a FSRQ at a redshift oi z = 0.859, a nd is one of 
the most intense and va r iable of these sour ces IjPian et al.l 
I2OO6I : IVillata et aLll2007l : lOupta et al1l2008l . and references 
therein). It is th e brightest and most variable blazar at 
GeV energies fe.g.'Verc ellone et al.ll2008l : lRaiteri et al.|[200i: 
IVerce llonc ct al. 2010). The multi-flare variability and SED 
study of the source in December 2 007 a nd De cember 2009 is 
report ed by iDonnarumma et aU l|2009l ') and iPacciani et al.l 
l|20ld '). respectively. The flaring behaviour of the source 
across the whole el ectromagnetic spectrum is studied by 
Ijorstad et al.1 (|2010|). They conclude that the emergence of 
a superluminal knot from the core yields a series of opti- 
cal and high-energy outbursts, and that the millimetre- wave 
core lies at the end of the jets acceleration and coUimation 
zone. 
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Figure 3. Modelled SED curves of the blazars PKS 1510-089 
and 3C 454.3. The last three points in all of these SEDs are almost 
certainly dominated by contribution from the disk emission, so 
we exclude these points while modelling the SEDs of these two 
FSRQs. 



5 DISCUSSION 

5.1 Statistical particle acceleration and 
log-parabolic spectra 

iMassaro et al.l IJ2004 ) showed that when the acceleration ef- 
ficiency is inversely proportional to the accelerating par- 
ticle's energy itself, then the energy distribution function, 
i.e., spectrum, approaches a log-parabolic shape. So, the log- 
parabolic spectra are naturally produced when the probabil- 
ity of statistical acceleration is energy dependent. According 
to this model, the curvature of emitting electron population 
(r) is related to fractional acceleration gain (e) as r oc 1/e. 
Also, Ep oc e follows a negative trend between Ep and b 
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(see lTramacere et al.ll2009l 'l. where Ep is the peak energy of 
observed SEDs. 

So, if the above analysis is correct, then one should ex- 
pect a negative trend between Ep or vp and h, where I'p 
is the peak frequency of observed SEDs. We analysed the 
correlation between the two quantities using our data and 
found a significant negative correlation between vp and b 
with rp = —0.74 and p-value = 0.00014, where rp is the 
linear Pearson correlation coefficient and the p- value for the 
null hypothesis of no correlation corresponds to a signifi- 
cance level > 99.98% (Fig. 4). So our result confirms that 
the spectral curvature parameter decreases as Ep moves to- 
wards higher energies and means that the blazars' spectra 
in our cases are fit well by the log-parabolic model. 

On the other hand, such a connection of log-parabolic 
spectra with acceleration may be also understood in 
the framework provided by the Fokker- Planck equation 
with momentum d i ffusion term (IKardashev et al.l 1 1963 : 
iMassaro et"aLl l2006l : iTramacere et al.ll2009l ). Thev showed 
that the log-parabolic spectrum results from a Fokker- 
Planck equation with a momentum diffusion term and a 
mono-energetic or quasi-monoenergetic injection where the 
diffusion ter m acts to broaden the s hape of the peak of the 
distribution. IKardashev et al.l (| 19621 ) showed that the curva- 
ture term, r ex: l/{Dt) where D is t he Diffusion term. Thi s 
relation leads to the following trend l|Tramacere et al.ll2009l ) 



ln{Ep) = 2 Ini-yp) + — 

00 



(5) 



So the both the momentum-diffusion term, D, and the 
fractional acceleration gain term, e, can explain the anti- 
correlation between Ep and b. 

An interesting check if the log-parabolic curve is actu- 
ally related to the statistical acceleration is the existence of a 
linear relation betwe en the two spectral parameters a and b 
(jMassaro et al.ll2004l V So we checked the correlation among 
these two quantities and found that they are significantly 
correlated, with rp — 0.61 and p- value = 0.003 (see Fig. 5). 
So, we consider that in our case the log-parabola type spec- 
tra are very likely to be characterized by a full statistical 
acceleration mechanism working on the emitting electrons. 

5.2 The cause of blazar SED changes 

There are several possibilities that could be invoked 
to explain the changes in the synchrotron spectra. In the 
following sub-sections, we will discuss some of them in detail. 



^ 



Log i^p,.,i, (IIz) 

Figure 4. Curvature (b) vs peak frequency {i^p)- A clear relation 
can be seen between the two quantities which is confirmed by a 
correlation analysis. The open symbols in this and subsequent fig- 
ures represent BL Lacs while the filled symbols stand for FSRQs. 
The straight line is the best linear fit with a slope, m = —0.14 
and constant, c = 2.22 {y = mx + c). 




5.2.1 Evolution of the particle energy density distribution 



The energy loss of the emitting particles is known to fol- 
low ^ = — C7^ where 7 is the Lorentz factor and C 
is a constant for both th e synchrotron and IC emission 
iRvbicki fc Lightmanlll979l ). The evolution of the relativis- 
tic particles can be described by ene rgy-dependent Fokker- 
Planck equation (see for instance iKembhavi fc Narlikan 
1 19991 . p.52) : 



dt 



d_[ dj 
'dj\dt'^ 



(7,i) +Q{l,t) 



(6) 



Figure 5. The spectral parameters, a and b. A clear relation 
can be seen between the two quantities which is confirmed by a 
correlation analysis. The straight line is the best linear fit with 
m = 0.47 and c = 0.96. 



where n(7,t) is the electron density distribution (a function 
of energy and time) and Q(7,t) is the time-dependent injec- 
tion rate. The solutions to the equation above are typically 
rather conrplex and can be obtained in terms of Green's 
functions for the general case. 

If the changes in the blazar SEDs are due to a grad- 
ual change of the electron energy density distribution due 



S 0.5 



Spectral Energy Distribution variation in BL Lacs and FSRQs 9 



^ 0.25 

o 



<1 



-0.4 -0.2 0,2 0,.4 0,6 



-2 



A(5 



Figure 6. Changes in the peak intensity vs. change in peak fre- 
quency for all blazars (both BL Lacs and FSRQs together) . There 
is no significant correlation found between the two quantities. 



Figure 8. Changes in the R band magnitude vs. change in 
Doppler factor for all blazars that have multiple SEDs in our 
sample. There is a significant correlation between the two quan- 
tities, with rp = —0.70 and p-value = 0.01. The straight line is 
the best linear fit with m = —0.23 and c = 0.29. 
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Figure 7. Changes in the peak intensity vs. change in Doppler 
factor for all blazars that have multiple SEDs in our sample. There 
is a clear relation seen between the two quantities. The straight 
line is the best linear fit with m = 0.07 and c = —0.07. 



to the synchrotron and IC losses, with no other injections 
meanwhile, then one should see both a total synchrotron 
energy decrease and mostly a decrease in 75 of the electron 
energy distribution (and respectively — the frequency of the 
synchrotron peak emission). On the other hand, freshly in- 
jected electrons between the observational epochs (their en- 
ergy should be higher than the average as it is expected to 
decrease in time) should reflect in the opposite behaviour 
(higher synchrotron emission and a peak blue-shift). So, in 



general one should expect to see a relation between the peak 
intensity and the peak frequency changes. To search for such 
a relation we used correlation analysis statistics (see Fig. 6). 
The correlation statistics revealed that there is no signifi- 
cant correlation between the two quantities as rp — 0.54 
corresponding to a p-value = 0.084, where rp is the lin- 
ear Pearson correlation coefficient and p is a null hypothesis 
probability value. As there is no significant relation found 
between the two quantities, perhaps electron energy density 
evolution/electron injection is not the primary driver of the 
SED changes. 



5.2.2 Change of the Doppler boosting factor 

Another effect that can modify the blazar synchrotron emis- 
sion is the change of the Doppler boosting factor, presum- 
ably d ue to a chang e eithe r the bulk velocity or the viewing 
angle. iRaiteri et al.l ([20101) also found that only geometrical 
(Doppler factor) changes are capable of explaining SED vari- 
ations between two epochs for BL Lac. The changes of the 
synchrotron peak intensity vs. Doppler factor variations for 
our sample of blazars are shown in Fig. 7. We can see a clear 
relation among the two quantities and it is confirmed by val- 
ues of rp = 0.51 and p-value = 0.05. We also searched for a 
correlation between change in R band magnitude vs. change 
in respective Doppler factor values. This would be relevant 
if <R> is a somewhat more representative quantity for the 
synchrotron emission rather than peak intensity. Again, we 
found significant correlation among ARmag and AS (see Fig. 
8). So we consider that Doppler factor changes are a strong 
driver that can be responsible for the SED changes. 
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Figure 9. Changes in the R band magnitude vs. change in mag- 
netic field for all blazars that have multiple SEDs in our sample. 
There is no significant correlation between the two quantities with 
rp = 0.11 and p-value = 0.75. 



5.2.3 Change of the Magnetic Field 

Since the value of the magnetic field, B is related to the 
location of the dissipation region, if the two SEDs arise 
from periods of different activity it is quite possible that 
the emission was produced by blob that dissipates at dif- 
ferent location in the jet. In that case we expect to find 
different B values for the different SEDs of our sources. The 
change in B values can also lead to variations in the flux 
of these sources which will be reflected in the two differ- 
ent SEDs of the sources. So we search for a correlation be- 
tween the change in <R> and that in B (Fig. 9) but we 
do not find any significant correlation among them for all 
of out blazars. However, we notice an apparent significant 
correlation among these two parameters for the BL Lacs 
alone, which is confirmed by correlation analysis {rp = 0.96, 
p = 0.002 for BL Lacs). So we conclude that changes in mag- 
netic field strength may be responsible for the SED changes 
in the case of BL Lacs. 



6 CONCLUSIONS 

We have carried out the radio to optical through mm, sub- 
mm and IR SED studies of a sample of ten blazars in- 
cluding five BL Lacs and five FSRQs, eight of which are 
LSPs and the other two ISPs. We modelled the SEDs of 
blazars using synchrotron spectra with log-parabolic distri- 
butions. We found a significant negative correlation among 
'-'Peak and the curvature term r, which implies that the ac- 
celeration efficiency of emitting electrons is inversely propor- 
tional to energy itself; however this correlation can also be 
explained by the momentum diffusion term in the solution 
of Fokker-Planck equation. Also, a significant correlation be- 
tween the two spectral parameters a and 6, implies that the 



log-parabolic curve is likely to be related to the statistical 
acceleration of emitting electrons. 

Of course our modelling has significant limitations. We 
are assuming that only one-zone dominates the emission at 
any given time, and this is unlikely to be an excellent ap- 
proximation. Although we strived for data obtained simul- 
taneously, this was often impossible to obtain; therefore any 
variations within the "high" and "low" states have been av- 
eraged over time bins ranging from days to two months. This 
lack of simultaneity could have vitiated our results but only 
appears to have added scatter to all of the Figures 4-9. 

We considered each likely factor that could be responsi- 
ble for the changes in observed SEDs of blazars. If the elec- 
tron energy density evolution governs the SED changes then 
one should expect a correlation between change in peak in- 
tensity vs change in peak frequency. Since we do not observe 
any such correlation we consider that probably the evolution 
of electron energy density cannot be responsible for the ob- 
served SED changes. Also for our entire sample of blazars, 
changes in <Kmag> are not correlated with the respective 
changes in B, so the change in magnetic field strength is 
probably not responsible for the SED changes we saw. How- 
ever, for the BL Lacs alone the SED changes may be driven 
by changes in B. We find that the change in Doppler factor 
is significantly correlated with the change in peak intensity 
as well as with <Kmag>. So it is reasonable to suggest that 
the change in Doppler factor (either bulk velocity or viewing 
angle) is the primary driver that governs the SED changes 
in short term variability (STV) in blazars. 
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